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Abstract 

In the framework of the three-band model for CuOi plane in high- 
temperature superconductors the spectrum of the spin-polaron hole 
exitation is investigated. The problem is treated taking into account 
the coupling of a local polaron with the antiferromagnetic spin wave 
with Q = (it, 7r).This leads to the essential changes of the lowest po- 
laron band ei(k) and the strong redistribution of the bare electron 
filling. 

PACS number(s): 74.72.-h. 75.10-b 

In order to understand the nature of high-temperature superconductors it is 
important to describe properly the motion of a hole in the C11O2 plane FJ] . 
This motion takes place on the antiferromagnetic (AFM) spin background of 
copper spins and it must be treated as a correlated motion of a hole coupled to 
spin excitations (a spin polaron). Usually the spin polaron is studied within 
the frameworks of the t-J model Jl[ and the three-band Hubbard model 
[0, |], m . In the previous works the spin polaron problem was studied mainly 
in the approximation of a small radius polaron ( an analogous of a Zang -Rise 
polaron). In the present work, for the first time, the long range order of the 
AFM background is taken into account by introducing an additional new 
spin polaron of an infinite radius -a bound state of a charge exitation and 
a spin wave with Q = (ir, 7r).We shall show that the introduction of such a 
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polaron leads to the essential decrease of lowest band filling by bare holes. As 
a result the area of the Fermi surface strongly increases and its form becomes 
more complex. 

The effective Hamiltonian of the three-band Hubbard model of the Cu0 2 
plane has the following form in conventional notations [§, ^ [5]: 

H = f + h + J, f = r £ *r 2 4 +a2)CT2 c R+ai)CTi (1) 

R, ai ,a 2 , 

CT 1 ,CT 2 

h — —h ^2 c R+a,<7 C R+a+b,<7 > J = 7^ ^2 S R S R+g . 

R,a,b,cr R,g 

Here the Cu0 2 plane is described by the square sublattice with lattice 
constant g and two O sites per a Cu unit cell; R — the vectors of Cu sites, 
R + a — are four vectors of O sites nearest to Cu site R, a = ia^, ia^, 
a x = g{\ify) a y — g(®, !)• We assume g = 1. In (P we take the notations: 
b are the nearest neighbor (NN) vectors for the oxygen sublattice; g = 2a 
and d = 2b are the first- and second- NNs for Cu sublattice; the operators 
c+ and X a0 create a hole with spin S = \ and spin projection | [o — ±1) 
at the O and Cu sites respectively and X R 1<72 are the Hubbard projection 
operators which are convenient for excluding doubly occupied Cu sites. 

The first term T in ([J) describes the effective hole hopping with amplitude 
from O to O sites through the intervening Cu sites. The term J correspond to 
AFM interaction between Cu sites, h represents the direct 0-0 NN hopping. 

Let us discuss the hole excitations with spin S = \ and the spin projection 
|. We shall restrict ourselves to a finite number of site operators Arj, (1 < 
j < 12), in each unit cell R. 

In order to treat the hole excitations in the framework of the spin-polaron 
concept we introduce for each unit cell R six site operators which describe 
the local polaron of small radius. 

^R,o-,l(2) = C R+a x (a y ),a ' ^R,cr,3(4) = ° X! 7^ R C R+a ^ a ^ 7 

7=±1 

^R,<r,5(6) = ° 7 X B,+g x ,(gy) C n+ ax ,( ay ),~f > O = -O (2) 

7=±1 

A + - 1 Vn ikR .-1.. ■ 
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A^ a j are the Fourier transforms of A^ a j. 

Let us mind, that in our previous investigations ]7) this basis of local 
spin-polaron operators lead to the proper description of the experimentally 
observed important features of the Cu0 2 plane hole spectrum: extended 
saddle point and isotropic band bottom. 

In this paper we want to investigate the role of the delocalized spin po- 
larons which correspond to the coupling of the local polarons to AFM spin 
wave with momentum Q = (71", 7r), so called Q-polarons usually, if the spin 
subsystem is found in the state with the AFM long range order then the 
average value of the amplitude (Sq) of the spin wave with q = Q (the Q- 
wave) has the macroscopic large value and has the properties analogous to 
the amplitude of a Bose particle with the zero momentum in the superffuid 
Bose-gas. As a result for many problems this amplitude can be treated as a 
c— number M. Then the coupling of the Q-wave to local electron states does 
not represent new states but leads, as mentioned above, to the mixing of the 
states with the momenta k and k+Q. This treatment is usually based on 
the widely used Neel type state of the spin subsystem with two sublattices. 

But for the S=l/2 spin system the quantum fluctuations are very im- 
portant and they lead to the spherically symmetric homogeneous Neel state 
at any small but finite temperature. In this background even at T = the 
average value (Sq) = and the mentioned above simple approach for the 
hybridization of k and k+Q states fails. In the homogeneous Neel state only 
(SqSq) can be treated as a macroscopic value. Then the coupling of a local 
polaron states to Sq corresponds to a new delocalized states. In order to take 
these states into account we introduce the additional six operators based on 
the basis of (||): 

A t*j = E tQr^r,^, = ° E 7^ +Q , 7 ,iS Q ? , (3) 

7=±1 7 
Q ya = ei QR S 7 * = N -l £ e iQ(R+R l)x 7^ j = j + 6, I = (1 6) . 

Ri 

In order to determine the spin polaron spectrum £j(k) of 12 quasiparticle 
bands we use the two time retarded matrix Green's functions Gij(t,\s) for 
the operators A k(Ti : 

Gy(f,k) = (A Kl (t)\Kj(°)) = "<©(*) ({ A kM . ( 4 ) 
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We solve the system of the equations of motion for Gij(u), k) by using the 
standard Mori-Zwanzig projection technique and restricting ourselves to the 
above chosen basis of operators {A k (Ji j} (0,0). Then the Green's functions 
and the spectrum are determined from the equations: 

(u - DK- 1 ) G = K , det\Ke(k) - D\ = 0. (5) 

A,i(k) = ({B Kl , A+ }) , K i4 = ({A Kl , A+J) , B Kl = [A kii , H . 

(6) 

The matrix elements of D and K are expressed both through short-range 
spin-correlation functions of Cu subsystem and the long-range order corre- 
lation function (SqSq). The Cu spin subsystem is described by the S = | 
Heisenberg model at T = 0. For the value of the correlation functions we use 
the results of || where this model was treated in the framework of spheri- 
cally symmetrical Green's functions theory. It is important that we take into 
account that ( 6/S'qSq) is a macroscopic value and is equal to a square of 
effective sublattice magnetization: 

(S Q S Q )= lim |(S R S R+Rl )| = M 2 (7) 

Note that the dependence of the spin polaron excitations spectrum on 
the long-range correlation function (SqSq) appears only due to treatment 
of the Q-polaron states @. Below we take the following numerical values 
of the spin-correlation functions: (SRSR +g ) = —0.3521, (SrSr +( i) = 0.229, 
(S R SR + 2 g ) = 0.2, M 2 = 0.0914. The matrix elements were calculated in the 
low hole doping limit, n << 1, n -is the total number of oxygen holes per 
unit cell. The detailed expressions of matrices D and K will be published 
elsewhere. 

As a result the Green's functions have a form 

12 ^(0 (\A 

GiM^^TT^jty (8) 
i=i u - £ i\ k ) 

In particular the value of Zh(k) = Z^^(k) +Z$ 2 J\t) corresponds to the 
number of bare oxygen holes with the fixed spin a and the momentum k 
in the state |k, a) of the lowest quasiparticle band £i(k). Let us mind that 
residue Z^j^(k) satisfy the sum rule 2 s -^(ij)(k) = 1, i = 1, 2. This means 
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that in this model the Luttinger theorem is not fulfilled and the maximum 
number of holes per cell is equal four despite the presents of twelve bands. 

The results of the most interesting two lowest bands £i(k) and £2 00 
for realistic values of model parameters J\ = 0.2r, h = 0.4r (below all 
energetical parameters are expressed in units r) are presented in Fig. 1(a). In 
the same figure the spectrum of the lowest of six bands £.i(k), calculated in 
the approximation of six operators (2) is shown. In Fig. 1(b) the spectrum 
£i(k) is presented by the equal-energy lines £i(k) = const. Let us mind that 
the Q-polarons may lead to a rather complex form of the Fermi surface. This 
circumstance can give a nontrivial behaviour of the Hall effect on doping, 
and even cause the inversion of the Hall constant if the Fermi energy is close 
to ei(k) = -4.5. 

As it is seen from Fig. 1(a), the inclusion of the Q-wave qualitatively leads 
to the decoupling of the lowest band of the local small polaron excitations 
and £i(k) is close to £i(k). This means that the main features of the lowest 
band excitations previously calculated in the local polaron approximation [[7j 
are preserved. 

The importance of the treatment of the Q-wave polarons may be seen if we 
discuss the filling of the lowest band by bare holes. In Fig. 1(c) the filling of the 
£i(k) an d £i(k) are shown, i.e., the residues -Zh(k) and Z_ h (k) (the underlined 
values correspond to local polaron approximation). One can see that the 
introduction of the Q-polarons leads to the essential decrease of hole filling 
along the lines X-M and M-N. This redistribution of the bare hole spectrum 
weight explaines the results of photoemission experiments when the "flat band 
region" is observed along the direction X-T and is invisible along the line X- 
M. Fig. 1(d) demonstrates that the local polaron concept (six operators (2)) 
leads to the strong decrease of the filling under the Fermi surface. This 
means the violation of the Luttinger theorem approximately in four times 
(the analogous effect was found in |lC]).But the inclusion of a Q-polaron 
states (3) leads to the addition essential reduction of filling, approximately 
in 1.5 times. The maximum £i(k) -band filling is equal to n = 0.22, and the 
smallness of this value justifies our low density approximation. 

Formula (3) points out that the Q-polaron A^ a7 ^ contains a bare hole 
state c^+q o-. It means that the residual Zq of the corresponding Green's 
function Gcj(k) = Grj + Gg,s are responsible for "shadow band" effect [[11 



We are to note that Q-polaron scenario reproduces the essential decrease 
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of the lowest band width with the decrease of the AFM constant J. Usually 
such an effect is obtained only in the self-consistent Born approximation [1^ . 

In conclusion we want to mention that the described above properties of 
the spin polaron are mainly preserved if we suppose that the spin subsystem 
has no long range order, but the spin correlation length L is large. Then, in 
the matrix elements of D and K we must replace the long-range correlation 
function (SqSq) by the following expression: 

E (SQ+qSQ+q) 5q « 1/L, L » 1. (9) 

|q|<<5g 
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Figure caption 



Fig.l: (a) The spectra £i(k) and £2(k) calculated for the basis (2), (3), 
and the spectrum £i(k) calculated for the basis (0). The spectra are given 
along symmetrical lines V — X — M — N — V and X — N — X; T = (0, 0); 
X = (tt, 0), (0, tt), M = (tt, tt); N = (tt/2, tt/2); 

(b) Spectrum equal-energy lines £i(k) = const; 

(c) Zh(k) and ^(k) — the number of bare holes (the residues of corre- 
sponding Green's functions) in the quasiparticle excitations for the spectra 
£i(k) and £_i(k); Zq(k) is the residue of the lowest pole £i(k); for the Green's 
function Gn t n(uj, k) +Cr 12,12 (<^, k), which characterize the "shadow band" ef- 
fect; 

(d) the dependence of the number of holes n per unit cell on the value of 
the Fermi surface area S (Sbz ~ the area of the first Brillouine zone): thick 
line — for the spectrum £i(k); dashed line — for the spectrum £_i(k); solid 
straight line — the case of the noninteracting particle filling. 
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